Fractional windings of the spinor condensates on a ring 
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We study the uniform solutions to the one-dimensional spinor Bose-Einstein condensates on a 
ring. These states explicitly display the associated motion of the super-current and the spin ro- 
tation, which give rise to fractional winding numbers according to the various compositions of the 
hyperfine states. It simultaneously yields a fractional factor to the global phase due to the gauge- 
spin symmetry. Our method can be applied to explore the fractional vortices by identifying the ring 
as the boundary of two-dimensional spinor condensates. 
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I. INTRODUCTION 

The Bose-Einstein condensation (BEC) of spinor 
atoms have been realized in the optical traps which 
confine the atoms regardless of their spin hyperfine 
states [2 [5]. Due to the spin degrees of freedom, there 
exist a variety of complicated structures such as mag- 
netic crystallization, spin textures as well as fractional 
vortices are allowed [3l-[l0]. The direction of the atom 
spin can change dynamically due to collisions between 
the atoms. The mean-field Hamiltonian of the spinor 
BEC is invariant under combined global gauge transfor- 
mation and spin rotation [/(l)phaso x '5'0(3)spin[5!. This 
spin-gauge symmetry reveals an interplay between super- 
fluidity and magnetism. A spatial variation of phase is 
always associated an internal rotation in the spin space, 
giving rise to a global phase in the condensate. This 
phase can be fractional times of 2tt which leads to a frac- 
tional vortex fT^H^ . 

The symmetries of the Hamiltonian are usually not 
completely broken in the ground state. According to the 
homotopy theory, the broken symmetry in the ground 
state specifies which type of the topological excitations, 
such as vortex, monopole, skyrmion, et al, can be hosted 
in that phase. [HI HO]- The set of operations under which 
the order parameter remains invariant constitutes an 
isotropy group to characterize the residual symmetry of 
the phase. Its conjugate class provides a label to classify 
the individual states. 

In this paper, we address the classification of fractional 
windings of the spin for an arbitrary spin-i^ {F = 1, 2, 3 
as example) BEC on a uniform ring. Non-uniform so- 
lutions are not appropriate to the present situation [TTj. 
We start from the simple plane wave solutions in which 
only two of the hyperfine states are nonzero. We pay our 
attentions to the symmetry and windings of the spin and 
the global phase around a circle on the ring. It shows 
that the states can be specified by a pair of quantum 
numbers {p,q), with the first number p {x2tt) specifying 
the global phase change around the ring and the second 
number q (x27r) the spin rotation. In general, p and q 
are fractional numbers. The corresponding state is said 
to have a fractional winding number. Furthermore, if we 



identify the boundary of a two-dimensional (2D) spinor 
BEC to the ID ring with large perimeter, then the ring 
state reflects possible topological defects residing in the 
2D system. It may provide a way to discern or create 
fractional vortices in the spinor BEC. 



II. SPIN-1 CONDENSATE 

The mean-field Hamiltonian of the spinor BEC is di- 
vided into two parts, 

H = Ho + Hs, (1) 
where the noninteracting part of the Hamiltonian is 

Ho^ [ drJ2^l[-^y' + V{rM^ (2) 

with V{r) the external potential. The interaction part of 
the Hamiltonian Hg is dependent of the spin of the BEC. 
For = 1 it reads [TU] 

H,^^Jdr[co:n'':+ci:P:], (3) 

where f are the three spin matrices and h is the density 
operator. The interaction strength cq and ci are related 
to the s-wave scattering length of total spin F = and 
F = 2 channels, respectively. 

The dynamical motion of the condensate is governed 
by the Heisenberg equation, which is explicitly written 
as the coupled Gross-Pitaevskii equations (GPEs), 

n— — 1 

where m = 0, ±1 indicate the hyperfine states and the 
spin polarization F = X]m('/'m|f iV'm)- The wave function 
should satisfy the periodic boundary conditions ijjm{L) = 
i-'m (0) on the ring. The stationary equations are obtained 
by substituting tpm{x,t) — -0„i(a;)e^*^*/'' into Eqs.(|4]), 
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= 



[--dl + (CO + C2)(|^i|2 + IV^oH + (CO - C2)\lP-l\^]^l + caVo'Ci 



1 



5^ + (Co + C2)(|^lP + l^-iP) + col^oPlV'o + 2c2^SVi^- 



[-\dl + (CO + C2)(|^-iP + \M^) + (CO - C2)m'']i>-i + C2^oVl* 



r 



(5) 



We consider the unform plane wave solution with ■00 
as 




(6) 



where 9 = ^-KxjL with L the perimeter of the ring. A 
and D are constants which can be directly calculated 
from Eqs.([5|, 



^--(— + — ) 
2co L'^^2cq 2c2^' 



2co y 2co 2c2 ' 



(7) 



In the limit of L — >■ oo, the corresponding ground state 
.^ground ^ (-^^ 23)'^ becomes a polar state ■0^°^'"' (with 
A~D^ 1/^/2) which has a continuous 5*0(2) symmetry 
and a discrete C2 symmetry perpendicular to the 5*0(2) 
symmetrical axis. The underlying reason is that _L — > 00 
implies that the kinetic energy is ignorable. For the finite 
ring sample (L is finite), A ^ D indicates only a discrete 
C2Z symmetry in Q since a super-current is excited along 
the ring. It is called the O2 state [T^. One notes that the 



A 

M 



42+j)2 is not simply 



mass circulation v^'^^^^'^dx 
quantized. 

As the full symmetry group of the Hamiltonian is 
J7(l)phase X 'S'0(3)spin, the mass flow of one hyperfine 
state will simultaneously induce a rotation in the spin 
space. It equivalently yields a global flow of the con- 
densate. The order parameter is described by V = 
y/n( = y/ne'^'^U{a, f3,j)(o, where (j) is the global gauge 
and U (a, (3, 7) = e~'^^''°'e~^^^^ e~'^^''^ is the rotational ma- 
trix. ^0 is the reference order parameter. By make a 
mapping a = —cj) — —9/2 from physical space into the 
order parameter space, we decompose the state ^ as 
t/j = e*^/2[/(-6'/2,0,0)(v4,0,£i)^, where the ID path is 
parameterized by the angular coordinate 9. In the state 
([6]), 9 increases 27r as one goes one circle along the ring. 

Figure 1 display the density distribution in the 
spherical-harmonic representation as^l ITUj. 



*(r,s) 



E 



■0™(r)YFm(s). 



(8) 



(b) 



FIG. 1: The spherical harmonic representation of the (|, |] 
state of the spin-1 BEG. (a) L = 1/2. (b) L = 2. 



The parameter are chosen as ao = 114aB, ^2 = lOSas, 
/i = 52.45 . It is seen that the the state exhibits the 0*22 
symmetry. The d vector describes the spin rotation when 
one goes along the ring. In this case, vector d rotates an 
angle of —it around the ring. 



Since tp{L) 



i/'(0), we get the global phase 



i X 27r and the associated spin-rotation angle a 
— ^ X 2n around the z-axis. Consequently, the state (6]) 
can be represented by a pair of quantum numbers ( 5 , f) ■ 
In general, an arbitrary current-carrying state on the ring 
can be specified by a pair of quantum numbers {p, q) with 
the p specifying the change of the global phase and q 
the change of the spin-rotating angle TU'. We observe 
that p and q are related to the mass current t;""^'*'* = 
^iV<jy-\f\Va) by p = ^ §\7cb-de and q ^ ^ §Va-de. 
Moreover, p and q are usually fractional numbers and the 
state is called to have fractional windings. Since the ID 
ring (with L — > 00) can be identified to the boundary 
of the 2D system, the state ip reflects a fractional vortex 
with winding numbers (|, |) residing in the 2D spinor 
BEC, due to the topological invariance. The vortex core 
may be flUed with normal phases |12j. 



3 



III. SPIN-2 CONDENSATE 

In this section we address the classification of fractional 
winding states of = 2 BEC on a ring. The interaction 
part of the Hamiltonian is [2^ 

Hs^^ dr[co : : +ci : : +C2iSo^o], (9) 



^■00 = 



where F+ = F* = 2(V'^Vi + V'-iV'-2) + y6(V'tV'o + 
^o*^_i), = 2(|V2p - \^-2?) + Hi? - and 



where Aqo = ■^{2ip2'4'^2 — '2'4'ii'-i+i^'a)- The stationary 
GPEs are 



(10) 



We consider the following forms of uniform solution in 
which only two of the hyperfine states are nonzero, 



V6 



C2 



[^2^^ ^ con± ciF2]'0±i + ci( — Fpi/)o + ^±^'±2) - -^AooV'j^i, 

1 /6 

+ Con + ^Ci(F+Vi + J^-^-i) + ^^00^0- 
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It is straightforward to calculate the constants A and D 
in each solution by substituting it into Eqs.(lO). The 



symmetry of the order parameters are analyzed by the 
same way as in the F — 1 BEC in the previous section. 



The ground state of Vi 



ground 



(A, 0,0, £1,0)^ is com- 



pared with cyclic state (^,0,0,^,0)-^ which has the 
tetrahedral symmetry T[9 . However, the symmetry of 
the Tpi state degenerates into the C32 symmetry by the 
same reason as in the spin-1 BEC that a super-current is 
excited. As 6 increases from to 2tt around the ring, the 
spin rotates an angle of — 27r/3 associated to a global 
phase of 27r/3 to yield i^iiL) = e^^/^^e^^-^Tr/a^^^O). 
Hence the ipi state has the (|, |) fractional winding num- 
bers with the symmetry, as shown in Fig. 2(a). 

The ground state of V'f """^ = (A, 0, 0, 0, D)^ IS com- 
pared to the biaxial nematic state (;^, 0, 0, 0, -^)'^ which 
has the I?4 symmetry (9j. Similarly, the symmetry of ■02 
has only the C4Z symmetry. As 6 increase from to 2tt, 
ij}2{L) = Q^'^ Q^h-^ I'^^i ^ Hence the state 1/12 has fractional 
winding numbers (^, ^), as is shown in Fig. 2(b). 



By the same way, we find that the ^3 state has the C22 
symmetry and can be specified by the winding numbers of 



round 



(A,0,Z?,0,-^)J 



IS 



{\-,\)- The ground state of 0| 
compared to the cyclic state (i, 0, 0, ■ However, 
the 04 state has only the Ciz symmetry due to the super- 
current in the ring. It is specified by the winding numbers 
of (0,i). 

As to the argument in the F = 1 BEC, from the states 
^1, '02, V'3, and -04 we conclude that the 2D F — 2 BEC 
may host (|, |), (|, |), (^, ^), and (0, ^) fractional vor- 
tices, respectivelv|10|. 



IV. SPIN-3 CONDENSATE 



The interaction part of the F — 3 Hamiltonian isjOlllOj. 



Hs = ^ dr[co : : +ci : P : -|-C2iJo^o+C3 ^1 ^2M^2m], 



(12) 
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FIG. 2: The order parameter of the spin-2 BEG for states specified as (a) ipi — >■ (|, |), (b) ?/'2 — > (§, j), (c) i/'s ~^ (!> §)' ^'^'^ 
(d) v^4^(0,i). 



where Aqo = ^(2'03?/;_3 
^(^V'±3V't2 - V5tp±2ip^i + ^fV'±i'0o), and i2±2 



7y(y f V'isV'Ti - ^Jf^±2^o + V2^±i)- The relevant 
physical parameters can be found in Ref. [S] . The station- 
ary GPEs are written as, 



'±1 



1 /6 c 



-^[5^2TlV43 - Vl5A2±l^;i - V20A2±2^o] 



2V2T 

[-191 + Con 



+ 



^;^[^^^2T2^;^a - \/5A2=plV'^2 - V3A2oi'^i + 2\/2A2±2V'±l + ^^2±l'0o]; 

\/3ci[i^_i/'-i +-P'+V'i] - ;^C2Aoo?/'o 
'*o + V2^2iV'i* + V2^2-iV':li - V20A22V^2* - ^20^22^12]- 



Ca 
2721 



[4^20^-0 



We consider the following combinations of uniform so- 
lution, 



(13) 
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(14) 
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FIG. 3: The order parameter of the spin-3 BEG for states 
specified as (a) Vi (|, j), (b) V2 ^ (|, |), (c) ^/'s ^ (5, |), 
(d) ^4 ^ (I i), and (e) ^5 ^ (i |). 



The constants A and I? in each state can be calculated 



directly by substituting the state into the GPEs (13). 



The symmetry of the order parameters can be ana- 
lyzed as in the previous sections. The phase diagram of 
the spin-3 BEG has been investigated in Ref. 9J. The 
state "01 has the C^z symmetry. As 6 increases from 
to 2tt around the ring, the spin rotates an angle of 
— 7r/4 associated to a global phase of 7r/4 which yields 

■>pi{L) = e''/^'e'f-''^^iJi{0). Hence the ipi state has the 
C4Z symmetry with winding numbers, as shown in 

Fig. 3 (a). The state -02 has the discrete C^z symmetry. 



When one goes around the ring, the spin rotates an an- 
gle of —2tt/5 associated to a global phase of An/ 5 which 
yields ^p2{L) = e'^^/5*e*-^-2^/5?Ai(0). Hence the state 02 is 
specified by a pair of fractional winding numbers ( | , | ) . 
Similarly, the state ips has the Cqz symmetry and is spec- 
ified by the winding numbers of (i, |). The state ■04 has 
the C^z symmetry and is specified by (|, |). The state 
^5 has the C4Z symmetry and is specified by ( 5 , | ) . 



■ii. ^ 



FIG. 4: The order parameter of the spin-3 BEG for states 
specified as (a) 06 -5> (0, |), (b) t/j7 -5> (0, |), (c) 08 -5> (0, i), 
and (d) V9 ^ (0, I). 

There is another class of combinations of the uniform 
state in which 0o 7^ 0. They are explicitly as follows, 
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(15) 



The symmetry and winding numbers are analyzed by the 
same way. The state tpe has the D^z symmetry and is 
specified by (0, |). The state ip7 has the Csz symmetry 
and is also specified by (0, |). The state ips has the D2 
symmetry and is specified by (0, i). Finally, the state 
"09 has the C2Z symmetry and is specified by (0, ^). This 
classification of states can be applied to the classification 
of fractional vortices in the 2D spin-3 BEG. From the 
states 01 ~ 09, we conclude that the 2D F ^ 3 BEG 
may topologically host the g, |, |, |, and ^ fractional 
vortices. 

V. SUMMARY 



We have studied the uniform solutions to arbitrary 
spin-F BEG on a ring. These states consist of various 



combinations of the hyperfine states. The symmetry of 
the states are analyzed. We demonstrated that the states 
can be specified by a pair of quantum numbers which de- 
scribe the global phase and the windings of the spin. Our 
method can help us to discern or create the fractional 
vortices in the 2D spinor BEG. 
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